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ABSTRACT 

V578 Mon is a system of two early B-type stars in the Rosette Nebula star-forming 
region (NGC 2244), and is one of only nine eclipsing binaries with component masses 
greater than 10 Mq whose physical parameters have been determined with an accu- 
racy of better than 3%. It is therefore a benchmark system for evolutionary and stellar 
structure models of newly formed massive stars. Combining our multi-band light curves 
spanning 40 yr with previous light curve data from the literature, we fit a model light 
curve that for the first time includes the effects of apsidal motion of the system. We mea- 
sure an apsidal period of 33.481q 06 y r - ^ s a consequence of incorporating the apsidal 
motion into the modeling of the system's orbital parameters, we determine an updated 
eccentricity of e = 0.07755+^^7 > 

which differs significantly from the value previously 
reported in the literature. Evidently, the inclusion of apsidal motion in the light curve 
modeling significantly affects the eccentricity determination. Incorporating these key 
parameters into a comprehensive model of the system's physical parameters — including 
internal structure constraints — will bring V578 Mon to the next level of benchmark 
precision and utility. 

Subject headings: binaries: close - binaries: eclipsing - individual: (V578 Mon) - stars: 
massive - stars: early type 



1. Introduction 

The analysis of apsidal motion in eccentric binary stars has been used for years to test stellar 
structure models. The periastron advance of an eccentric binary system is a direct consequence 
of the finite size of the stellar components, and of tidal interactions be tween them. Consequently, 



apsidal motion serves as a measure of the internal structure of stars (jSternd Il939l ). Specifically 
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measurement of apsidal motion in eccentric bin ary systems allow stringe nt tests of the internal 
structure constant hi predicted by theory (e.g. I Claret fc Gimenea 12010). Apsidal m otion also 
provides a test of General Relativity outside of our Solar System (e.g. IWolf et al.ll2010l ). 



The photometric variability of the 2.408 day period, eccentric, massive eclipsing binary (EB) 
V578 Mon (HDE 259135, BD+4°1299), co mprising a Bl type primary star and a B2 type secondary 
star, was first identified in the study by iHeiserl (|l977l ) of NGC 2244 within the Rosette Nebula 
(NGC 2237, NGC 2246). The absolute dimensions of V578 Mon have been deter mined from three 



season s of Stromgren ubvy photometry and one season of radial- velocity data by iHensberge et al 



(|2000l . hereafter H2000). An analysis o f the m etallicity and evolutionary status of V578 Mon was 
undertaken by iPavlovski &: Hensbergd (|2005l ) and H2000. The masses and radii of V578 Mon 
determined from these data are 14.54 ± 0.08 M and 10.29 ± 0.06 Mq, and 5.23 ± 0.06 R and 
4.32 ± 0.07 Rq, for the primary and secondary respectively (H2000). 

These masses and radii are accurate to better than 3%, making V578 Mon one of only nine 



EBs w ith Mi > M2 > 10 Mq and with sufficient accuracy to be included in the I Torres et al 



(|2010l ) compilation of benchmark-grade EBs. However, of these, V578 Mon is the only eccentric 
EB lacking an apsidal motion measurement. 

The classical theory of tides and General Relativity both predict that a close, eccentric system 
such as V578 Mon will experience a certain amount of periastron advance. The angle of periastron, 
u, is given by the equation u(t) = wo + tbt, where ojq is the angle of periastron at the reference 
epoch HJDo, and t is the time since HJDo- The apsidal period is given by U = 360°/cO x P, 
where P is the orbita l period and Cj is the apsidal motion in deg cycle -1 . The AAVSO research 
note by IHeiserl (|2010l ) presented the first long term photometry of V578 Mon that was used to 
identify apsidal motion in the system and to update its ephemeris. The author estimated a value of 
U ~ 30.4 yr from 14 independent primary eclipse minima measured over a 40 year period. Here we 
present an ana lysis of the apsidal mo tion of V578 Mon using the state-of-the-art EB modeling soft- 
ware PHOEBE (IPrsa &: Zwitterl 120051) . which updates and extends the venerable Wilson-Devinney 



code rtWilson &, Devinney 



1971; [Wilson 



mine d via eclipse timings (e.g. iGimenez &: Quintana 



19791) . Traditionally, apsidal motion of EBs has been deter 



1992 



Gimenez fc Bastero 



1995 



Wolf fc Zeida 



20051 ). However, the apsidal motion of an eccentric EB causes not only the eclipse timings to vary; 
it also causes the shapes and depths of the eclipses — as well as the morphology of the out-of-eclipse 
portions of the light curve — to vary over time. Therefore, in principle a full light curve model takes 
into account more of the apsidal information encoded in the light curve data, and should yield an 
extremely precise measure of the apsidal motion. 

Using our own light curve data spanning 40 yr together with previous light curves from the 
literature, we measure an apsidal period for V578 Mon of 33.48^ 'q 6 yr. Furthermore, as a conse- 
quence of including the apsidal motion for the first time into the analysis of the orbit of V578 Mon, 



we report an updated eccentricity of e = 0.07755 



+0.00022 

-0.00027 



, which differs from the previous literature 



value of 0.0867±0.0006 (H2000). These fundamental orbital parameters set the stage for follow-up 
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analyses to determine the internal structure of the stars in V578 Mon for the first time, and to 
re-determine the stellar radii, for detailed tests of stellar evolution models with this benchmark 
system. 

In £}2j we present the photometry used in this paper. In Sj3j we perform the light curve analysis 
of all photometry. In SJU we present the apsidal period and orbital eccentricity along with an error 
analysis. We conclude in ^5] with a discussion and summary. 



2. Data 



The available time-series photometry of V578 Mon covers nearly 40 yr and more than one full 
apsidal motion period. A summary of the various light curv e epochs, incl uding filters and observing 
facilities used, is presented in Table [TJ Photometry from iHeiserl (|2010l ) includes multiband light 
curves spanning 1967-2006 from the 16-in telescope at Kitt Peak National Observatory (KPNO) 
and from the Tennessee State University(TSU) -Vanderbilt 16-in Automatic Photoelectric Telescope 
(APT) at Fairborn Observatory. The KPNO Johnson UBV light curves comprise 725 data points 
spanning 1967-1984 with average formal uncertainties per data point of 0.004 mag. The APT 
Johnson BV light curves span 1994-2006 and consist o f 1783 data points with formal uncertainties 
per data point of 0.001 mag for B and 0.002 mag for V (|Heiserll2010l ). Light curves from H2000 span 
1991-1994 from the 0.5-m Stromgren Automatic Telescope (SAT) at La Silla, with 248 data points 
in each of the ubvy filters and average formal uncertainty per data point of 0.003 mag (H2000). 
Table Q] lists these formal uncertainties, <ro, as reported by the original authors. However, from our 
light curve fits (see below) we found that these formal errors were in most cases underestimated. 
Thus we also report as a in Table Q] the uncertainties that we ultimately adopted for each light 
curve (see §3.21 for details). 



3. Light Curve Analysis 

Here we determine oj with the light curve modeling program phoebe using all of the available 
photometry (Table [1]). We perform the light curve fitting in two separate approaches. In the first 
approach (b |3. 1 j) , we determined Co by finding the linear change in uj(t) from fits to individual light 
curve epochs. In the second approach ( M3.2j) . we determine Co by finding a global light curve solution 
simultaneously to all of the light curve data in Table Q3 

All fixed parameters are listed in Table [2l For fixed parameters in both approaches, we adopt 
the spectral type of the primary to be B1V, which implies T\ = 30000K (H2000, and references 
therein). We adopt gravity brightening (gi, c^) and surface albedos (A\, A%) to be 1, as appropriate 
for stars with radiative envelopes. The rotational synchronicity parameters {F\, F2) are fixed at 
1.13 ± 0.03 and 1.11 ± 0.03 for the primary and secondary, respectively, based on the vs'mi and 
radii determined by H2000. Our limb darkening coefficients follow the square-root law for fully 
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radiative stars (|Claretll2000l ). In all cases we also adopt the semi-major axis (a) and the mass ratio 
(q = M2/M1) from H2000 because these parameters are determined principally from the radial- 
velocity curves. We do not include star spots in any of the light curve modeling; as discussed below, 
the variations in the light curve data are fully reproduced through the effects of apsidal motion. 



3.1. Fits to individual light curve epochs 

To obtain an initial, simple estimate of the apsidal period, we first determine uj by finding the 
linear change in w from the individual Johnson BV light curve epochs 1973-1976, 1994-1995, and 
2005-2006 (Tabled]). We chose these epochs because they span one full apsidal period, and because 
these light curves were obtained using the same instrument and filter set. Johnson BV light curve 
epochs 1973-1976 are a portion of the KPNO 1967-1984 light curve listed in Table U 

The free parameters are the angle of periastron w, the inclination of the orbit i, the secondary 
temperature T%, and the surface potentials and O2. The starting values for these parameters 
are from H2000 (see Table [2]). Note that we fit for the surface potential of the stars (i.e., Q oc 
i? -1 ) as well as T2, in order to give the light curve model full freedom to fit the eclipse widths 
and depths, however we regard the best-fit values of these parameters as preliminary. We fixed 
e = 0.077551q'qqq27, as determined in the global solutio n to al l light curve data described in §3.21 



We set the orbital period , P, t o the value from iHeiserl (|2010l ). and we used the primary minima 



eclipse times from IHeiserl (|2010l ) as the HJDq values for each light curve epoch. 



The resulting best-fit uj at specified HJDo's and formal errors from each light curve epoch 
are given in Table EJ A least-squares line fit to these uj vs. orbital cycles initially gave a reduced 
chi square x^ c d = X 2 /N = 19.0, where N is the number of data points. This indicates that the 
uncertainties on the individual uj values are underestimated. This is not surprising because the 
formal uncertainties on uj here do not take correlations with other parameters into account, as we 
do in our global light curve solution below. Therefore we scaled the uncertainties in Table [3] to 
achieve Xred = 1- The slope of the fitted line gives uj = 0.0718 ± 0.0012 deg cycle -1 and thus an 
apsidal period U = 33.06 ± 0.58 yr. 



3.2. Global fit to all light curve data 



Next we perform a global light curve solution fit simultaneously to all epochs and filters of 
available light curve data (see Table [TJ). The free parameters are e, ujo, i, T2, Qi, and uj. The 
initial values for these parameters are from H2000 except for uj, for which we used uj = 0.0780 
deg cycle -1 calculated from the est i mated apsidal period from IHeiserl (|2010l ). We again fixed P 
and HJDq to the values from IHeiserl (|2010l ). 



We perform two global light curve fitting iterations. For the first iteration we used the formal 
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photometric uncertainties, do in Table HJ which yielded a total X^. ed = 3.25, indicating underesti- 
mated photometric errors. Thus for the second iteration we scaled the photometric errors to make 
the x 2 °f each light curve to equal the number of data points N of the light curve. The scaled 
photometric errors a listed in Table [1] reflect more realistic, conservative values given the quality 
of our photometry. Effectively this causes X, r ed ®^ each light curve to be approximately unity. This 
yielded a total x 2 e d = 0-97 for the final fit. 

The best fit global solution is plotted over multiple light curve epochs in Figures [21 SI 
andO The resulting apsidal motion is u> = 0.07089^q'qqq^3 deg cycle -1 , giving an apsidal period 
U = 33.48^QQg yr. The quoted uncertainties are from a detailed \ 2 analysis (see £ j4.3p . In Table [2j 
we list the full set of system parameters resulting from the global light curve model fit. 

Note that the changes in the shapes, depths, and timings of the primary and secondary 
eclipses — as well as the changes in the out-of-eclipse portions of the light curves — are due to apsidal 
motion effects and are very well reproduced by the model. These variations over time due to apsidal 
motion are clearly demonstrated in Figure [61 which displays the global light curve solution for the 
Johnson B light curve epochs 2005-2006, 1999-2000, 1995-1996, 1994-1995, and 1973-1976. 



4. Results 

4.1. Apsidal period of V578 Mon 

We have determined the apsidal motion, 6j, of V578 Mon via model fits to the available light 
curve data. Calculating the linear change in uo(t) by fitting individual light curve epochs spaced over 
one apsidal period results in to = 0.0718 ±0.0012 deg cycle" 1 and U = 33.06 ±0.58 yr. Calculating 
d) via a global light curve solution results in to = 0.07089tQQQQ^3 deg cycle -1 and U = 33.48^0'og yr. 
The two values are consistent within the uncertainties. We adopt the latter value as it incorporates 
the full dataset, is more precise, and furthermore our error analysis accounts for correlations in the 
fitted parameters (see below). 

We prefer the approach of fitting the full light curves because it incorporates all available light 
curve data in Table HJ not just the eclipse timings. As shown across multiple observing seasons in 
Figures El [3l and El the shapes, depths, and timings of the primary and secondary eclipses change 
with time due to uj. Furthermore, the light curves show variation in the out-of-eclipse data over 
time which are very well reproduced by the light curve model. The model light curves do not 
include star spots, showing that the out-of-eclipse variations in the light curve of V578 Mon are 
indeed a manifestation of apsidal motion. 



- 6 - 



4.2. Orbital eccentricity of V578 Mon 

The eccentricity from our light curve model fit incorporating apsidal motion is e = 0.07755^0 ooos 
(Table[2]), which differs significantly from the previously reported value of 0.0867± 0.0006 (H2000). 

To further investigate the H2000 eccentricity, we perform a light curve fit within PHOEBE using 
only the SAT photometry (the same data used in H2000). We set Cj to zero and a, q, wq, i, T 2 , f^i 
and to values from H2000. The only free parameter is e. This phoebe fit converges to a light 
curve solution with e = 0.0867, reproducing the e found by H2000. Evidently, accounting for the 
effects of apsidal motion yields a different e. The e, Qi and O2 are correlated parameters in light 
curve analysis, meaning that our significantly different e could yield different radii for the stars. 
The new, tentative radii we compute are approximately 5.14 Rq and 4.70 Rq for the primary and 
secondary as compared to the literature values of 5.23 ±0.06 Rq and 4.32 ±0.07 Rq (H2000). Thus 
there is an indication that the updated e together with the Cj newly reported here may result in a 
significantly different R2. This will be the subject of an in-depth analysis in a forthcoming paper. 



4.3. Uncertainties on Cj and e 

In order to determine realistic uncertainties on Cj and e, we performed a detailed analysis of 
the x 2 space around the best fit values. We varied Cj over 0.0702-0.0716 deg cycle" 1 with step 
length 0.00006 deg cycle -1 , and we varied e over 0.0765-0.0787 with step length 0.000084. For each 
of the 625 combinations of Cj and e, we recomputed the global light curve fit as before. 

Figure [7] shows the resulting contour plot of \ 2 for e vs. Cj. Contours are drawn at A% 2 
values corr esponding t o la, 2a, and 3a uncertainty for a A% 2 distribution of two parameters 



of interest (jPressI 119881 ). The contour shapes indicate that Cj and e are not strongly correlated 
given our analysis. Moreover, Figure [7] clearly demonstrates that the previously reported value of 
e = 0.0867 ± 0.0006 (H2000) lies well beyond the 3a contour and can therefore be ruled out with 
very high statistical significance. 



5. Discussion and Summary 



The accurate apsidal period of 33.481q 06 y r an< ^ updated orbital eccentricity of e = 0.07755to'ooo27 
of V578 Mon underscores the value of a long time baseline of photo metric observations for ecc entric 
eclipsing binary stars. The traditional eclipse timing method (e.g. iGimenez fc Basterdll995l ) uses 
the timing of the primary and secondary eclipse of an eccentric binary star system to calculate e 
and Cj among other parameters. In this paper, we use the eclipse timings but also the changing 
shapes, widths, and depths of the primary and secondary eclipse due to Cj, as well as the light curve 
variations in the out-of-eclipse phases to determine realistic constraints on Cj and e. Furthermore, 
we demonstrate that including the apsidal motion parameter in light curve fitting can affect the 
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eccentricity measurement of eclipsing binary systems. 

The apsidal motion of V578 Mon can provide an accurate test o f theor etical calculations of 



the internal structure constant (A^)- For example, IClaret &: Gimenea ([20101 ) find generally good 
agreement between the theoretically predicted and measured hi when they consider EBs with radii 
to ±2% accuracy. Comparing against EBs with very accurately measured radii is critical, because 
the theoretically predicted hi is highly dependent upon the stellar radii (&2 oc R 5 ). Currently, 
only 18 eccentric EBs have stellar radii and Cj measured with sufficient accuracy to stringently test 
theoretical internal structure models. Of these, only EM Car, V478 Cyg, V453 Cyg, and CW Cep 
have masses greater than 10 Mq, well constrained lj, and radii accurate to ±2% ([Claret &: Gimenez 



2010J, and references therein). A detailed reanalysis of V578 Mon will yield accurate radii and, 



combined with our u = 0.07089^q'qqq^3 deg cycle 1 , will yield an accurate &2 with which to test 
theoretical internal structure models of high mass stars. 

There is still much to know about high mass stellar evolution, especially at young ages, mak- 
ing V578 Mon an important testbed for stellar evolution models given the large amount of precise 
photometry and radial velocity data on the system. To date, there are only ni ne EBs with compo - 



nent masses greater than 10 Mq whose masses and radii are accurate to ±3% ([Torres et alj|2010l ). 
Stellar evolution models for stars with masses greater than 10 Mq thus remain poorly constrained 
by EBs. Future reanalysis of V578 Mon will include a preci s e cali bration of high-mass evolution 



models similar to the work on V453 Cyg by lSouthworth et al.l ([2004]). Finally, V578 Mon's location 



allows for precise age and distance constraints of the Rosette Nebula, similar to previous work on 
V578 Mon by H2000. A follow-up paper will incorporate the accurate orbital parameters newly 
determined here in order to re-determine all of the physical properties — including internal structure 
parameters — of this important, benchmark EB. 

We acknowledge the support of the NSF REU program in Physics & Astronomy at Vanderbilt 
University, and NSF grants AST-0849736 and AST-1009810. 
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Table 1. V578 Mon Light Curves 



Observatory 


Year 


Filter 


O"0 


a 


N 








(mag) 


(mag) 




1 KPNO 


1967-84 


Johnson U 


0.004 


0.016 


251 






Johnson B 


0.004 


0.012 


256 






Johnson V 


0.004 


0.013 


217 


2 SAT 


1991-94 


Stromgren u 


0.0029 


0.0067 


248 






Stromgren b 


0.0023 


0.0046 


248 






Stromgren v 


0.0023 


0.0054 


248 






Stromgren y 


0.0030 


0.0053 


248 


3 APT 


1994-95 


Johnson V 


0.0037 


0.0022 


260 






Johnson B 


0.001 


0.0040 


254 


APT 


1995-96 


Johnson V 


0.002 


0.0035 


95 






Johnson B 


0.001 


0.0037 


96 


APT 


1999-2000 


Johnson V 


0.002 


0.0058 


259 






Johnson B 


0.001 


0.0078 


246 


APT 


2005-06 


Johnson V 


0.002 


0.0036 


284 






Johnson B 


0.001 


0.0044 


283 





Note. — He-inch telescope at Kitt Peak (KPNO) 
2 0.5 m telescope at La Silla (SAT) 

3 TSU-Vanderbilt 16-inch telescope at Fairborn University (APT) 
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Table 2. V578 Mon Parameters 





Free Parameters 


H2000 Value 


This study 


Apsidal Motion, Cj [deg cycle -1 ] 




0-070891™ 


Apsidal Period, U [yr] 




33-48±S:S8 


Angle of Periastron, wq [deg] 


153.3 ±0.3 


159.8 ±0.33 


Eccentricity, e 


0.0867 ± 0.0006 


nTT^^+O- 00018 

U.Uf 'OO_0.00026 


Inclination, i [deg] 


72.58 ±0.30 


71.67 * 


Primary Surface Potential, f^i 


5.02 ± .05 


5.15 * 


Secondary Surface Potential, 


4.87 ± .06 


4.56 * 


Secondary Temperature, T2 [K] 


26400 ± 400 


26100 * 




Fixed Parameters 


Ephemeris, HJDo [d] 




2449360.6250 


Orbital Period, P [d] 




2.4084822 


Primary Surface Temperature, T\ [K] 




30000 ± 500 


Systemic Velocity, 7 [kms -1 ] 




34.9 ±0.1 


Semi-major Axis, a [R@] 




22.03 ±0.04 


Mass Ratio, q = M 2 /M\ 




0.7078 ± .0002 


Primary Synchronicity Parameter, F\ 




1.13 ±0.03 


Secondary Synchronicity Parameter, F2 




1.11 ±0.03 





Note. — * These parameters are preliminary. 



Table 3. u(t) values from fitting individual light curve epochs 



Light Curve 




HJD 


u 






[days] 


[deg] 


Johnson BV 2005- 


-06 


2453741.73648 


104.00 ±0.26 


Johnson BV 1994- 


-95 


2449738.75624 


348.10 ±0.69 


Johnson BV 1973- 


-77 


2443105.90927 


140.8 ± 1.6 
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Fig. 1. — Variation of angle of periastron oj determined by model fits to individual light curve 
epochs. The circles represent u) values as determined from light curves at the epochs 1974-1977, 
1994-1995, and 2005-2006 (Table [3]). Uncertainties on the individual uj are smaller than the 
plotting symbols. The dashed line represents a linear fit using the equation oj(t) = ujq + Cot, where 
t is time since the reference epoch, HJDo- This results in to = 0.0718 ± 0.0012 deg cycle" 1 , and 
U =33.06 ±0.58 yr. 
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Fig. 2.— Representative fits to light curves from 2005-2006, 1999-2000, 1995-1996 and 1994-1995 
in the Johnson B passband from global fits to all light curve data, offset for clarity (see §3.2p . The 
residuals to the fits (O — C) are shown above. 
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Fig. 3. — Same as Fig. [21 but showing Johnson V band light curves and fits. 
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Fig. 4. — Same as Fig. [21 but showing Stromgren ubvy light curves and fits. 
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Fig. 5. — Same as Fig. [21 but showing 1973-1977 Johnson UBV light curves and fits. 
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Fig. 6. — Representative light curve model fits in the Johnson B passband from our global light 
curve solution (see §3.2n . Note the variation of the shapes, depths, and timing of the primary and 
secondary eclipses. The out-of-eclipse portions of the light curves also vary, due to the apsidal 
motion, not star spots. Solid lines represent fits to the 2005-2006 and 1973-1976 epochs. The 
dashed, dash-dot, and dotted lines represent fits to the 1999-2000, 1995-1996, and 1994-1995 
epochs, respectively. 
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Fig. 7. — Constraints on the apsidal motion Co and eccentricity e described in §4.31 The contours 
represent A% 2 = X 2 ~Xmin v& l ues of 2.30, 6.170, and 11.80 (corresponding to 1, 2, and 3 a confidence 
intervals). The black dot represents our best fit solution. The resulting best fit and uncertainties 
are u = 0.07089±g;gggf| deg cycle" 1 and e = 0.07755±g;ggg|| . 



